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1.   Introduction. 

r  25 1 

Deterministic  Lanchester-type  equations  of  warfare  (see  TAYLOR  and  BROWN   '  , 

[27L 
WEISS    )  play  an  important  role  in  military  operations  research  for  developing  insights 

[41 
into  the  dynamics  of  combat  (see,  for  example,  BONDER  and  FARRELL     ,  or  BONDER  and 

HONIG   ) ,  even  though  combat  between  two  opposing  military  forces  is  a  far  more  complex 

[9] 
random  process.   The  classic  Lanchester  theory  of  combat  (see  DOLANSKY   )  considered 

constant  attrition-rate  coefficients.   New  operations  research  techniques  for  fore- 
casting temporal  variations  in  fire  effectiveness  (caused  by,  for  example,  changes  in 
force  separation,  combatant  postures,  target  acquisition  rates,  firing  rates,  etc.) 
have  generated  interest  in  variable-coefficient  combat  formulations.   Unfortunately, 
the  resultant  differential  equations  are  not  well  studied. 

In  this  paper  we  present  a  mathematical  theory  for  predicting  battle  outcome 
from  initial  conditions  without  explicitly  computing  force-level  trajectories  for  variable 
coefficient  Lanchester-type  equations  of  modern  warfare  for  combat  between  two  homogeneous 
forces  (see  Note  1) .   The  determination  of  conditions  on  initial  values  that  predict 
(in  the  sense  of  necessary  and/or  sufficient  conditions)  force  annihilation  in  such 
Lanchester-type  combat  leads  to  some  new  mathematical  problems  in  the  theory  of  ordinary 
differential  equations.   This  force  annihilation  problem  may  be  viewed  as  either  a 
problem  of  determining  the  asymptotic  behavior  of  the  solution  (depending  on  given 

initial  conditions)  or  a  problem  of  determining  the  range  of  the  quotient  of  two  linearly 

[25] 
independent  solutions  to,  for  example,  the  X   force-level  equation     (see  Note  2). 

In  either  case,  the  classic  ordinary  differential  equation  theories  (see ,  for  example, 
HILLE1    ,  INCE*-    ,  OLVER^   ')  are  inadequate  to  supply  all  the  answers  sought.   We 
show  that  questions  of  force  annihilation  can  be  reduced  to  the  study  of  certain  "expo- 
nential-like" Lanchester  functions  and  may  be  simply  answered  by  examining  certain 
inequalities  involving  the  initial  conditions  and  by  possibly  consulting  tabulations 
of  new  special  functions  that  are  suggested  here.   Our  general  results  apply  to  a  wide 
class  of  attrition-rate  coefficients  (namely,  those  that  yield  continuous  force-level 
trajectories) . 


Thus,  in  this  paper  we  provide  a  general  theoretical  framework  for  determining 
force  annihilation  without  explicitly  computing  force-level  trajectories  for  variable- 
coefficient  Lanchester-type  equations  of  modern  warfare.   [Other  modes  of  battle  ter- 
mination are  briefly  discussed.]   We  introduce  a  canonical  auxiliary  parity-condition 
problem  for  such  determinations.   New  exponential-like  general  Lanchester  functions 
arise  from  the  solution  to  this  problem,  and  tabulations  of  these  would  facilitate 
force-annihilation  prediction.   Different  mathematical  forms  for  attrition-rate  coeffi- 
cients lead  to  different  auxiliary  parity-condition  problems.   Our  theory  is  applied 
to  general  power  attrition-rate  coefficients:   exact  force-annihilation  results  are 
obtained  for  cases  of  "no  offset"  (modelling,  for  example,  weapon  systems  with  the 
same  maximum  effective  range);  and  although  qualitative  results  are  obtained,  future 
computational  work  is  required  for  quantitative  results  for  cases  of  "offset"  (modelling 
for  example,  weapon  systems  with  different  maximum  effective  ranges) . 

2.   Lanchester 's  Classic  Formulation. 

ri3i 

F.  W.  LANCHESTER     hypothesized  in  1914  that  combat  between  two  military  forces 
could  be  modelled  by 

dx/dt  =  -ay,  dy/dt  =  -bx,  (1) 

with  initial  conditions 

x(t=0)  =  xQ,  y(t=0)  =  yQ,  (2) 

where   t  =  0  denotes  the  time  at  which  the  battle  begins,   x(t)   and  y(t)   denote 
the  numbers  of  X  and  Y  at  time   t,   and  a  and  b  are  nonnegative  constants  which 
are  today  called  Lanchester  attrition-rate  coefficients  and  represent  each  side's  fire 
effectiveness.   Lanchester  (see  McCLOSKEY     for  his  influence  on  operations  research) 
considered  this  model  (1)  in  order  to  provide  insight  into  the  dynamics  of  combat  under 
"modern  conditions"  and  justify  the  principle  of  concentration  (see  Note  3) .   We  will 
accordingly  refer  to  (1)  as  Lanchester 's  equations  of  modern  warfare.   Various  sets  of 

physical  circumstances  have  been  hypothesized  to  yield  them:   for  example,  (a)  both 

[27] 
sides  use  aimed  fire  and  target  acquisition  times  are  constant  (see  Weiss    ) >  or  (b) 


both  sides  use  area  fire  and  a  constant  density  defense  (see  BRACKNEY   ) .   Other  forms 

of  Lanchester-type  equations  appear  in  the  literature,  but  we  will  not  consider  these 

T91  T221 

here  (see  Dolansky1  J  and  TAYLOR1   J). 

From  (1)  Lanchester  deduced  his  classic  square  law 

b(X2-x2(t))  =  a(y2-y2(t)).  (3) 

Consider  now  a  battle  terminated  by  either  force  level  reaching  a  given  "breakpoint" 

op 

(see  Note  4):   for  example,   Y  wins  when  xf  =  x(tf)  =  Xgp  =  f^   xQ  but  yf  >  yBp 

BP 
=  f  yA,   where   t,.,   xc ,   y  c     denote  final  values  and  x^..   denotes  X's  breakpoint 
I   (J  r     I     I  or 

BP 
which  is  given  fraction  f    of  his  initial  strength.   It  follows  from  (3)  that 

A 

Y  wins  «  xQ  <  [{l-(f^P)2}/{l-(f|P)2}]1/2  /aTb  yQ,  (4) 

which  for  a  f ight-to-the-f inish  becomes  the  classic  result 

Y  wins  f ight-to-the-f inish  »  x.  <  /a/b  y  .  (5) 

Since,  unfortunately,  no  relationship  similar  to  (3)  holds  in  general  for  variable 
attrition-rate  coefficients,  we  observe  that  (4)  may  also  be  obtained  from  the  time 
history  of  the  X  force  level 

x(t)  =  {(x0-y0v^7b)exp(/ab  t)  +  (xQ+yQ/a/b)exp(-/ab  t)}  /2  ,  (6) 

BP 
via  determining  the  time  for  X   to  reach  his  breakpoint  (i.e.   x(t=tv  )  =  xDT)) 

t|P  =  (l//ab)£n({-xBp+/x^p+y^a/b-x^}/{y0^7b  -  xQ})  ,  (7) 

BP     BP 
and  requiring   t_  <  tL.  .   The  key  result  for  obtaining  (7)  is  that  one  of  the  two 

linearly  independence  solutions  to  the  X  force-level  equation  d2x/dt2  -  abx  =  0   is 

the  reciprocal  of  the  other.   For  a  f ight-to-the-f inish  (7)  becomes 

tJJ  =  {l/(2v^b)}Jln({y0^7b  +  x0}/{y0^7b-x0}),  (8) 

where  t   denotes  the  time  to  annihilate  the  X  force.  We  observe  that  (5)  is  an 

A 

immediate  consequence  of  (8) . 


In  many  applications  (see  Section  3  below) ,  one  is  interested  in  whether  the 

battle  will  be  terminated  within  a  given  time   t  .   In  this  case  x„  <  /a/by    is  a 

g  0      J0 

necessary  condition  for  X  to  be  annihilated  and  annihilation  occurs  when   t  £  t  . 

X   g 

Thus,  determination  of  whether  force  annihilation  will  occur  within  a  given  time  involve 
consulting  a  tabulation  of  a  transcendental  function,  here  the  natural  logarithm  (see 
also  (10)  below).   Similar  results  hold  for  other  fixed  force-level  breakpoints. 
The  time  history  of  the  X  force  level  may  also  be  written  as 

x(t)  =  xn  cosh  v^ab  t  -  y  /a/b  sinh  /ab  t .  (9) 

[251 
Taylor  and  Brown     take  (9)  as  their  point  of  departure  for  a  mathematical  theory  for 

solving  variable-coefficient  formulations.   (7)  does  not  follow  directly  from  (9),  but 


(5)  does  via  x(t=t  )  =  0  and 


-1. 


t*  =  (l/v^b)tanh  1(/b7a"x()/y0),  (10) 

since  the  range  of  the  hyperbolic  tangent  is   [0,1]   for  nonnegative  arguments. 

The  purpose  of  this  paper  is  to  generalize  the  above  to  the  general  case  of 
variable  attrition-rate  coefficients.   We  use  (6)  as  our  point  of  departure  and  base 
our  development  on  the  observation  that  (5)  follows  directly  from  (6) ,  since  the  second 
term  in  brackets  is  always  positive  and  goes  to  zero  as   t -*■+«>.   [We  will  ignore  the 
physical  impossibility  of  negative  force  levels  in  developing  results  like  (5).]   Thus, 


by  (6)   lim  x(t)  =  -<»  «  x_  <  /a/b  yn ,   whence  follows  (5). 
t->  °° 

3.   Variable  Attrition-Rate  Coefficients. 

The  pioneering  work  of  S.  Bonder        on  methodology  for  the  evaluation  of 

military  systems  (in  particular,  mobile  systems  such  as  tanks)  has  generated  interest 

in  variable-coefficient  Lanchester-type  equations  and  has  led  to  improved  operations 

\2]    [3] 
research  techniques  for  the  prediction  of  such  coefficients  (see  BONDER       ;  back- 

[251 
ground  and  further  references  are  given  in  Taylor  and  Brown    ).   Thus,  we  consider 

dx/dt  =  -a(t)y,  dy/dt  =  -b(t)x,  (11) 

where  a(t)   and  b(t)  denote  time-dependent  Lanchester  attrition-rate  coefficients. 


These  coefficients  depend  on  such  variables  as  firing  doctrine,  firing  rate,  rate  of 
target  acquisition,  force  separation,  tactical  posture  of  targets,  etc.  (see  reference 
4).   Without  loss  of  generality,  we  may  take  a(t)  ■  k  g(t)   and  b(t)  =  k,  h(t)  ,  where 
g(t)   and  h(t)   denote  time-varying  factors  such  that  a(t)/b(t)  =  k  /lc  =  constant 
for   g(t)  =  h(t).   We  will  also  refer  to  (11)  as  the  equations  for  a  square-law  attri- 
tions process ,  since  an  "instantaneous"  square  law  holds  even  when  a(t)/b(t)   is  not 

f  261 
constant  (see  TAYLOR  and  PARRY1   J;  also  references  21,  22,  and  24). 

A  large  class  of  combat  situations  of  interest  can  be  modelled  with  the  following 

attrition-rate  coefficients  (see  reference  4) 

a(t)  =  ka(t+C)y  and  b(t)  =  kb(t+C+A)V,         (12) 

where  A,C  ^  0.   We  will  refer  to  these  coefficients  as  general  power  attrition-rate 

[25] 
coefficients .   The  modelling  roles  of  A  and   C  are  discussed  in  Taylor  and  Brown. 

We  will  refer  to  C  as  the  starting  parameter,  since  it  allows  us  to  model  (with 
y,v  ^  0)   battles  which  begin  within  the  maximum  effective  ranges  of  the  two  systems. 
We  will  refer  to  A  as  the  offset  parameter,  since  it  allows  us  to  model  (again,  with 
y,v  £  0)   battles  between  weapon  systems  with  different  effective  ranges.   For  example, 
let  us  consider  Bonder's    constant-speed  attack  on  a  static  defensive  position  (see 
also  references  22  and  25) .   Then  we  have 

dx/dt  =  -a(r)y  -  -aQ(l-r/Ra) yy ,  dy/dt  =  -B(r)x  =  -B0(l-r/R  ) Vx,    (13) 

where   u,v  ^  0  and  R   denotes  the  maximum  effective  range  of  Y's  weapon  system 
(i.e.   a(r)  =  0  for  r  >  R  ),  and  these  parameters  are  given  by  (we  assume  that 


R6  *  Ra  *  V 


A  =  (R  -R  )/v,  and  C  =  (R  -Rn)/v,  (14) 

pa  a   u 


where  R   denotes  the  battle's  opening  range  and  v  >  0  denotes  the  constant  attack 
speed.   By  considering  (14)  and  Figure  1,  the  reader  should  have  no  trouble  in  under- 
standing our  terminology  for  A  and   C.   In  this  model  (13)   u ,   for  example,  is  used 
to  model  the  range  dependence  of  Y's  attrition-rate  coefficient  (see  Figure  2). 
Negative  values  of   u  and  v  are  discussed  in  reference  25.   Range  is  related  to  time 

by  r(t)  =  R  -vt.   Since  zero  force  separation  is  reached  at  time   t  =  R^/v*   the 
0  r  g    0 
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above  model  (13)  provides  interest  in  predicting  battle  termination  within  a  given 
time   t   (see  Section  2  above) . 

Almost  all  previous  work  on  the  variable-coefficient  equations  (11)  has  developed 
infinite  series  solutions  for  force-level  trajectories  or  represented  these  by  tabulated 
functions  (see  reference  25,  in  particular  Section  3).   Relatively  little  attention  has 
been  given  to  determining  the  qualitative  behavior  of  solutions  to  (11)  (such  as  predic- 
tion of  battle  outcome)  without  explicitly  computing  battle  trajectories.   Bonder  and 

[4] 
Farrell    (see  Note  5),  however,  have  considered  force  annihilation  within  a  given 

time.   Using  comparison  techniques  from  the  theory  of  ordinary  differential  equations 

tot 
(see,  for  example,  CODDINGTON  and  LEVINSON1   ) ,  they  obtained  a  rather  strong  sufficient 

condition  for  the  special  case  of  the  model  (13)  with  u  =  v  =  1  and  A  >  0. 

4.   A  Mathematical  Theory  for  Conditions  of  Force  Annihilation. 

Motivated  by  the  constant-coefficient  results  (see  Note  6) ,  we  introduce  the 
exponential-like  general  Lanchester  functions  E  ,  E  ,  EL,,   and  ELT,   defined  by 

dE+/dt  =  v^Tir"  a(t)E^  with  E^(t=tQ)  =  1/Q, 


+  .  .   _+. 


dE^dt  =  /k7k^b(t)Ex  with   E^(t=tQ)  =  1 


(15) 


and 


dE, 


7dt  =  S\l\   a(t)Ey  with  E~(t>t0)  =  1, 


d^/dt   =   S\J\  Mt)Ex  with     Ev(t=t0)    =  Q; 


(16) 


X  Y  X 

where   t_  =  max(t  ,t  ) ,   and   t   denotes  the  largest  finite  singular  point  (see  Taylor 

and  Brown    ;  also  p.  69  of  Ince    )  (see  Note  7)  on  the  t-axis  for  the  X  force- 
level  equation 

d2x/dt2 -{[l/a(t)]da/dt}dx/dt-a(t)b(t)x  =  0.  (17) 

By  assumption,  then,   a(t)   and  b(t)   are  positive  continuous  functions  V  t  >  t,,. 

Both  Ex(t)  and  Ex^t^  satisfy  (17).  Since  they  are  a  fundamental  system  of 
solutions,  we  may  use  them  to  construct  all  solutions  to  (17).  Thus,  the  solution  to 
to  (17)  which  satisfies  the  initial  conditions  (2)  is  given  by 
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X(t)  =  {[x()EY(t=0)-yir7k^y0Ex(t=0)]E^(t) 

+  [x0E^(t=0)  +  ^7k^"y0Ex(t=0)]Ex(t)}/2,     (18) 
and  from  (11)  we  obtain 


y(t)  =  {[y0Ex(t=0)-»1cb/ka  x()EY(t=0)]EY(t) 

+  [y04(t=0)  +  v/V^IxoEY(t=0)]EY(t:)}/2'         (19) 

where  we  have  made  use  of  the  easily  verifiable  fact  that  (see  reference  25) 

Ex(t)E^(t)  +  Ex(t)E^(t)  =  2  Vt.  (20) 

Without  further  restrictions  on  a(t)   and  b(t),   the  systems  (15)  and  (16) 
might  not  have  a  solution:   it  might  be  that  one  or  more  of  the  general  Lanchester 
functions  is  unbounded  at   t_.   (see  Section  4  of  reference  25).   For  example,  unbounded- 

ness  occurs  for  power  attrition-rate  coefficients  with  u  +  v+2  =  0  and  A  =  0   (see 

[251 
equation  (37)  of  Taylor  and  Brown    ).   Thus,  we  introduce 


CONDITION  (A) : 


t*t0. 


a(s)ds   and 


ft 

b(s)ds  are  bounded  for  all  finite 


'o  fco 


Then,  we  have 

THEOREM  1:   Condition  (A)  is  a  necessary  and  sufficient  condition  for  (15) 
and  (16)  to  have  a  continuous  solution  for  all  finite   t£t  . 

PROOF :   Sufficiency  follows  by  Theorem  6.4.2  on  p.  226  of  Hille^   J  (see  also  p.  64 

of  LEE  and  MARKUS    ^).   To  prove  necessity,  it  suffices  to  consider  E  (t)  .   If  a 

a 

continuous  solution  exists,  then  from  (15)  we  have  dE  /dt £  /lc/k  a(t),  whence 

a(s)ds,   and  the  theorem  follows.   Q.E.D.   Let  us  henceforth 
t, 


E*(t)  £l/Q+/ku/k 
X  b  a 

0 

assume  that  Condition  (A)  holds 


Let  us  now  consider  how  to  choose  the  general  Lanchester  functions  E   and  E 

A  A 

so  that  they  play  the  roles  of  a  decaying  exponential  and  an  increasing  exponential, 

9 


respectively,  in  the  solution  (18)  to  the  variable-coefficient  X  force-level  equation 

(17)  [whence  our  notation  of  E   and  E  ].  We  recall  the  constant-coefficient  result 

X         X 

(6)  and  its  consequence  (5)  obtained  using   lim  exp(-/ab  t)  =  0  and  exp(-/ab  t)  >  0. 

t-*-H» 

For  any  positive  value  of  Q  it  is  clear  that  E  (t)  ,  which  satisfies  (15),  grows 

x 

without  bound  just  as  an  increasing  exponential  does.   We  will  now  show  how  to  choose 
E   so  that  it  corresponds  to  a  decaying  exponential.   Similar  statements  hold  for 

A. 

Ey  and  E„.   Considering  (16),  we  see  that  we  should  choose  E   and  E   to  remain 
positive  for  all  t  so  that  by  (16)  they  continuously  decrease.   Furthermore,  we  will 
be  able  to  specify  such  behavior  for  E   and  E   by  our  selection  of  the  parameter 

A  I 

Q  in  the  initial  conditions  for  (16) . 

The  solution  E  (t)  ,  E  (t)   to  (16)   depends  continuously  on  the  parameter  Q 

A.        X 

of  the  initial  conditions  (see,  for  example,  Hille    ).   We  denote  this  dependence 
by  E~(t;Q),  E^(t;Q).   Let  Q*  =  Q*(a(t) ,b( t) )   denote  the  value  of  Q  such  that 

E~(t;Q=Q*),  E^(t;Q=Q*)  >  0  for  all  finite  t  ^  tQ.  (21 

It  follows  from  (16)  that 

lim  E~(t;Q*)  =  lim  E^(t;Q*)  =  0.  (22 

It  is  intuitively  obvious  that  such  a  Q*  exists,  and  we  will  prove  its  existence 

in  some  particular  cases  below.   As  we  shall  see,  knowledge  of  Q*  provides  valuable 

information  about  the  qualitative  behavior  of  force-level  trajectories  for  the  Lan- 

chester-type  equations  (11) .   Let  us  refer  to  the  problem  of  determining  Q*  such 

that  (22)  holds  as  the  auxiliary  parity-condition  problem.   Unless  explicitly  stated 

otherwise,  for  convenience  we  will  denote,  for  example,   E  (t;Q*)   as  E  (t) . 

X  X 

Comment  1:   For  a  constant  ratio  of  attrition-rate  coefficients,  i.e. 

a(t)  =  kah(t),  and  b(t)  =  kfah(t) ,  (23 

where  h(t)  denotes  the  common  time-varying  factor  of  the  two  coefficients,  it  readily 

follows  from  the  results  given  in  references  4  and  20  that  Q*  =  1  and 

10 


E*(t)  =  E^(t)  =  exp{ip(t)},       and       Ex(t)  =  E^(t)  =  exp{-^(t)},     (24) 


where  ip(t)  =  /k  le     h(s)ds 


[251 
Comment  2:   By  Theorem  1  of  Taylor  and  Brown    ,  the  solution  (18)  simplifies  to  the 

form  of  (6)  only  if  (23)  holds. 

The  above  exponential- like  general  Lanchester  functions  may  be  related  to  Taylor 

[251 
and  Brown's     hyperbolic-like  general  Lanchester  functions.   Let 

Cx(t)=X;L(t),    Sx(t)=x2(t),    CY(t)=y1(t),    and   Sy(t)=y2(t),     (25) 

where  x1 ,  x~  ,  y1 ,   and  y„  denote  the  general  Lanchester  functions  introduced  by 

Taylor  and  Brown.   Then,  similar  to  the  well-known  relationships  between  the  hyperbolic 

and  exponential  functions,  we  have 

Cx(t)  =  {Q*Ex(t)+Ex(t)}/2,  Sx(t)  =  {Q*Ex(t)  -Ex(t)}/(2Q*), 

+  +  (26) 

C^t)  =  {Q*EY(t)+Ey(t)}/(2Q*),  SY(t)  =  {Q*Ey(t)-EY(t)}/2. 

The  determination  of  Q*  will  be  slightly  simplified  for  general  power  attrition- 
rate  coefficients  (12)  by  considering  a  modified  auxiliary  parity-condition  problem. 
For  this  purpose  we  introduce  the  new  independent  variable 

rt 


=  K/lc  /k      a(o)da, 


(27) 


fc0 


and  define  s  =  s(t=0)^0   for   t  £0.   K  is  an,  at  present,  undetermined  parameter. 
It  will  be  chosen  so  that  a  more  convenient  canonical  system  of  differential  equations 
arises  in  the  modified  auxiliary  parity-condition  problem.   By  Condition  (A) ,  the 
transformation  (27)  is  well  defined  for   t^tn-   It  has  an  inverse   t(s)  ,   since  a(t)  > 
V  t  >  t„.   Letting 

ex(s)  =  KEx(t(s)),   e+(s)  =  E+(t(s)),   e~(s)  =  E~(t(s)),   and   e~(s)  =  E~(t(s))/K, 
the  substitution  (27)  transforms  (15)  through  (17)  into 

de*/ds  =  e*  with   e*(s=0)  =  1/Z, 

+  +  +  (29) 

de^/ds  =  Ks)e  with   e^(s=0)  =  1, 

11 


de„/ds  =  -£>  with   e  (s=0)  =  1, 

de^/ds  =  -Ks)ex  with   ^(s=°)  ■  z> 


(30; 


and 


where  for  any  Q 


and 


d2x/ds2-I(s)x  =  0,  (31 


Z  =  Q/K,  (32] 


I(s)  =  ({b(t)/kb}/{a(t)/ka>)/K2,  (33] 

[12] 
is  the  invariant  of  the  normal  form  (31)  (see  p.  119  of  KAMKE    )  and   t  =  t(s)   by 

(27).   The  parameter  K  will  be  chosen  to  simplify  the  form  of  I(s).   In  our  later 

work  the  equation  (31)  will  be  easier  to  analyze  than  (17) . 

We  will  refer  to  the  problem  of  determining   Z*  =  Z*(a( t) ,b(t) )   such  that 

e~(s;Z=Z*),  e~(s;Z=Z*)  >  0  for  all  finite  s  £  0,  (34] 

A  I 

as  the  modified  auxiliary  parity-condition  problem.   By  (32)  we  then  have  that  Q*  =  KZ; 

We  also  observe  that  the  solution  to  (31)  which  satisfies  (2)  may  be  expressed 
in  terms  of  these  exponential-like  general  Lanchester  functions  for  s  ^  s-   as 

x(s)  =  {  [x0Ke^(s=s0)  -  J\f\   y0e~(s=s0)]e^(s)/K 

+   [x0Ke^(s=s0)  +  \/\  y0e^(s=s0)]e^(s)/K}/2,  (35] 

and   similarly 


y(s)    =   {[yQex(s=s0)  -  *\/*T  xQKeY(s=s0)]eY(s) 

+    [y0ex(s=s0)  +  /k^7F"x0KeY(s=s0)]eY(s)}/2,  (36; 

where  by    (20)  and    (28) 


e*(s)e~(s)  +ex(s)e^(s)    =2  Vs 


(37; 


From  our  choice  of  Q*  such  that  (21)  and  (22)  hold,  we  can  immediately  infer 

the  behavior  of  the  solution  (18)  to  (11)  as   t  ->-  +°°  and  similarly  for  y(t)  (19). 

Thus ,  we  have 

11 


THEOREM  2:  lim     x(t)    =  -«     if   and   only   if 

x0E~(t=0;Q*)    <    AiJ\  yQE~(t=0;Q*). 

Equivalently,  we  may  state 

THEOREM  3:   Consider  combat  between  two  homogeneous  forces  described  by 
(11).   Assume  that  (11)  applies  for  all  time  and  that  Y   "wins"  when 
x(tf)  =  0  with  y(tf)  >  0.   Then,  Y  wins  if  and  only  if  x  E~(t=0;Q*)  < 
^7k^y0E-(t=0;Q*). 

Thus,  we  see  that  tabulations  of  such  new  exponential-like  general  Lanchester  functions 
E  (t;Q*)  and  Ey.(t;Q*)  would  facilitate  force-annihilation  prediction.  By  our  choice 
of   t_   only  one  such  tabulation  is  necessary  for  given  attrition-rate  coefficients 

a(t)   and  b(t)   (see  Note  7).   Alternatively,  we  may  express  the  force-annihilation 

[251 
condition  of  Theorem  3  in  terms  of  Taylor  and  Brown's     hyperbolic-like  general 

Lanchester  functions  (see  equations  (25)  and  (26)).  We  have  then 

THEOREM  3':   Consider  combat  between  two  homogeneous  forces  described  by 
(11).   Assume  that  (11)  applies  for  all  time  and  that  Y   "wins"  when 
x(t  )  =  0  with  y(t  )  >  0.   Then,   Y  wins  if  and  only  if  x  {Q*Cy(t=0)  - 
Sy(t=0)}  <  \/\   y0(Cx(t=0)  -Q*Sx(t=0)}. 

As  an  immediate  corollary  to  Theorem  3  we  have 

COROLLARY  3.1:   For   t  =  0,  Y  wins  if  and  only  if  x  Q*  <  A   /k^  yQ . 

One  shortcoming  of  our  above  development  is  that  Theorem  3  and  its  corollaries 
are  basically  existence  theorems  for  the  annihilation  of  one  side  by  the  other  at  some 
unknown  future  point  in  time.   As  the  enemy  initial  force  level  decreases  towards 
parity  (i.e.  equality  holding  in  (38)),  the  time  required  to  annihilate  a  force  becomes 

larger  and  larger.   There  is,  in  fact,  no  limit  to  how  large  it  may  become.   Moreover, 

[4] 
there  is  a  large  class  of  tactically  significant  battles  (see  Bonder's    constant- 
speed  attack  on  a  static  defensive  position  discussed  in  Section  3  above)  which  has  a 

13 


built-in  time  limit,  denoted  as   t  .   Hence,  it  would  be  desirable  to  have  a  method  for 

g 

determining  from  initial  conditions  (without  explicitly  computing  the  entire  force-level 
trajectories)  whether  or  not  force  annihilation  will  occur  within  a  given  finite  time 
t  .   For  our  general  model  (11),  Theorem  3  tells  us  that 


x0F^(t=0;Q*)  <  \/\   yQE~  (t=0;Q*),  (38) 

3.  3.  3. 

is  a  necessary  condition  for  x(t„)  =  0  with   t„  ^  t  ,  where  t,,  denotes  the  time  at 

X  X    g  X 

which  the  X  force  is  annihilated.   Motivated  by  the  well-known  constant-coefficient 
result  (8) ,  one  intuitively  sees  that  determining  whether  or  not   t  £  t   will  require 
the  appropriate  tabulations  of  new  transcendents  (i.e.  new  functions).   Different  such 
functions  arise  from  different  fundamental  systems  chosen  to  construct  the  solution  to, 
for  example,  the  X  force-level  equation  (17).   Let  us  now  investigate  which  fundamenta 
system  of  solutions  is  the  most  useful. 

We  begin  our  investigation  by  outlining  for  the  new  exponential-like  general 
Lanchester  functions  how  to  determine  from  initial  conditions  (without  explicitly  com- 
puting the  force-level  trajectories)  whether  or  not  force  annihilation  will  occur  within 
a  given  finite  time  t  .   Looking  at  (18)  and  setting  x(t)  =  0,  we  see  that  we  must 

O 

solve  for  E  (t)/E  (t)  =  n(t)   (see  Note  8).   For  any  other  fundamental  system  (i.e.  pai 

A        A 

of  solutions,  we  must  still  solve  for  such  a  quotient  (cf.  the  constant-coefficient 
results  (8)  and  (10)).   Thus,  our  force-annihilation  determination  requires  the  use  of 
tabulations  of  the  quotient  of  two  linearly  independent  solutions  to,  for  example,  the 
X  force-level  equation.   Our  question  is  now  which  quotient  is  the  most  useful  (i.e. 
which  fundamental  system  yields  the  most  useful  quotient)  .  We  will  show  that  the  quo- 
tient of  two  exponential-like  general  Lanchester  functions  is  not  numerically  satis- 
factory for  such  determinations,  whereas  the  quotient  of  two  hyperbolic-like  general 
Lanchester  functions  is. 

Motivated  by  the  result  for  a  constant  ratio  of  attrition-rate  coefficients  (23) 
that  n(t)  =  exp{2ij>(t)},  we  use  the  following  notation  for  the  n(t)   of  the  X  force- 
level  equation 
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E2X(t)  =  Ex(t)/EX(t)'  (39) 

with   E„  (t)   being  similarly  defined.   Assuming  that  (38)  holds,  we  see  from  (18) 

that  if  x(t=ty)  =  0,   then 

A 


E2x(tx>  =  {^J\  y0Ex(t=0;Q*)+xo4(t=0^)}/ 


{   a^b  y0Ex(t=0;Q*)  -x0EY"(t=0;Q*)}. 
We  will  prove  below  that  E_  (t)  =  E_  (t;Q*)   is  a  strictly  increasing  function  of   t 

with  initial  value  1/Q*  at   t  ■  tft.   Consequently,  the  inverse  function  E    (£)  = 

+ 
L_  (£)   is  a  well-defined  function  of   £  V  £  z   [l/Q*,+°°).   Hence,  we  may  write  that  the 

time  to  annihilate  X  is  given  by  (see  Note  9) 

\  =  E2x"1({/1V\  yoEX(t=0;^  +x04(t=0;Q*)}/ 

i\/\  y0E~(t=0;Q*)  -xQE^(t=0;Q*)}).  (40) 

Considering  the  above  developments,  one  may  show  that  except  when  (23)  holds,  it  is 
not  possible  to  determine  in  the  manner  of  (40)  the  battle-termination  time  when  X's 


breakpoint  is  positive,  i.e.   x^p  >  0, 


We  now  show  that  E.  (t)   (defined  by  (39))  is  a  strictly  increasing  function. 


We  readily  compute  using  (15),  (16),  and  (20)  that 


+ 


dE2x/dt  =  2/1^/^  a(t)/{Ex(t)}2     with  E2x(t=tQ)  =  1/Q*,  (41) 

whence  follows  the  monotonicity . 

[Sometimes  it  is  convenient  to  express  results  in  terms  of  the  modified  exponen- 

+       —       +  — 

tial-like  general  Lanchester  functions   ey(s) >   ev(s) »   ^(s) »   and   eyCs) »  which  are 

based  on  the  transformed  time  variable  s   defined  by  (27)  .   The  above  results  are 

readily  extended  to  this  case.   In  particular,   e   (s)  =  e  (s)/e  (s)   is  a  strictly 

increasing  function  of   s   for  s  ^  0,   e„  (s=0)  =  K/Q*,   and   lim  e~  (s)  =  +°°.] 

^X  _u_  *■*>. 

Thus,  (I)  determination  of  Q*  =  Q*(a(t) ,b(t) )  ,   and  (II)  tabulation  of 

E.  (t)   and  E~  (t)   would  allow  one  to  determine  whether  or  not  force  annihilation 

occurs  in  such  battles  with  finite  time  limit  without  explicitly  computing  the  entire 

15 


force-level  trajectories  (see  Figures  3,  4,  and  5  of  reference  25).   Unfortunately, 
there  is  a  serious  drawback  to  considering  E_  (t)   and  E_  (t)  :   accurate  tabulations 
are  difficult  (in  fact,  essentially  impossible  for  large  values  of   t)   to  generate, 
since  both  functions  are  basically  increasing  exponentials  so  that  any  error  in  their 
initial  value  1/Q*   (which,  in  general,  can  be  only  numerically  approximately  deter- 
mined) becomes  tremendously  magnified  over  time. 

We  may  develop  numerically  satisfactory  functions  for  prediction  of  force  annihi- 
lation within  a  given  finite  time,  however,  by  considering  the  hyperbolic-like  general 

r  25 1 

Lanchester  functions  of  Taylor  and  Brown    .   Let  us  therefore  define 

Tx(t)  =  Sx(t)/Cx(t).  (42) 

The  functions  T  (t)   and  Tv(t)   are  analogous  to  the  hyperbolic  tangent,  to  which 

they  reduce  for  a  constant  ratio  of  attrition-rate  coefficients.   Considering  (25)  and 

equation  (16)  of  reference  25,  we  see  that,  for  example,   Tv(t)   does  not  depend  on  Q*, 

X 

since  S  (t)   and  C  (t)   do  not.   Thus,   T  (t)   and  T  (t)   are  numerically  suitable 
XX  X  Y 

for  determining  whether  or  not  force  annihilation  will  occur  within  a  given  finite  time. 

[251 
Using  the  results  given  in  Table  I  of  Taylor  and  Brown    ,  we  readily  compute 

that 


dTx/dt  =  \l\   a(t)/{Cx(t)}2     with  Tx(t=tQ)  =  0.  (43) 

Hence,   T  (t)   is  a  strictly  increasing  function,  and  its  inverse  T    is  well  defined. 
A  X 

Let  us  now  establish  an  upper  bound  for  T  (t).   By  (26)  and  (42),  we  have 

A 

Tx(t)  =  (l/Q*){Q*Ex(t)-Ex(t)}/{Q*Ex(t)+Ex(t)}, 
whence  it  follows  that  for   t  ^  t_ 

0  *  T(t)  <  1/Q*,     with  lim  T(t)  =  1/Q*.  (44) 

X  X 

t-H-°° 

Thus,  our  current  investigation  has  yielded  important  information  about  the  asymptotic 
behavior  of  hyperbolic-like  general  Lanchester  functions . 
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3.  3. 

To  determine   t   such  that  x(t  )  =0,   we  write  the  solution  to  (17)  which 

A  A 

r  251 
satisfies  the  initial  conditions  (2)  as1   J   x(t)  =  x  {C^(t=0)C  ( t)  -  S  ( t=0)S  (t) } 


-  y  /kTk^  {C  (t=0)Sx(t)  -  Sx(t=0)Cx(t)}  and  find  that  when  (38)  holds,   t   is  given  by 
CX  =  TX1({x0CY(t=0)+y0,/V1S  Vt=0)}/ 


tyQAa\   Cx(t=0)+xQSY(t=0)}).  (45) 


We  observe  that  by  (26)  the  argument  of  the  inverse  function  T    in  (45)  belongs  to 

A 

the  range  of  T    (see  (44))  when  (38)  holds  (see  Theorem  3').   For  t_  =  0,   (45)  simpli- 

a    -1   / 

fies  to   t  =  T   (/k,/k  x-Vy-)  .  Whether  or  not  force  annihilation  occurs  within  a 

given  finite  time   t    then  depends  on  whether  or  not   t,,  £  t  .   Thus,  (I)  determination 

g  X    g 

of  Q*,   and  (II)  tabulation  of  Ty(t)   and  ^(t)   would  allow  one  to  determine  (with- 
out explicitly  computing  the  force-level  trajectories)  the  time  at  which  a  side  is 
annihilated. 

5.  Application  to  Power  Attrition-Rate  Coefficients. 

Let  us  now  apply  the  above  general  theory  to  (11)  with  the  general  power  attrition- 
rate  coefficients  (12).   We  observe  that  in  this  case   tn  =  -C,  where  C  ^  0.   In  order 
that  Condition  (A)  holds  we  must  have  u,v  >  -1. 

As  we  have  seen  in  Section  4,  our  theory  of  force-annihilation  prediction  depends 

on  knowing  Q*,   the  solution  to  the  auxiliary  parity-condition  problem.   For  the  power 

attrition-rate  coefficients  (12),  it  is  more  convenient,  however,  to  determine  Q*  via 

the  modified  auxiliary  parity-condition  problem  (30)  (see  also  (34)).   Hence,  we  apply 

the  transformations  (27)  and  (28)  to  (16).   For  the  coefficients  (12),  equations  (30) 

take  the  form 

dev/ds  =  -e  with  ev(s=0)  =  1, 

X        Y  X  (46) 

de^/ds  =  -sB(l-h-/sa)Vex     with  e^(s=0)  =  Z, 

where  the  parameter  K   in  (27)  is  given  by  K  =  ( /k  k  /(u+1))  P   ;   and  we  have 
p  =  (u+l)/(u+v+2),   q  =  1-p,  a  =  l/(u+l),   8  =  (v-y)/(u+l),   and 
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2/(u+v+2) 


Y  =  A(/k  tc  /(M+D)  /vm    '.   For   u,v  >  -1,   we  have  0  <  p,  q  <  1. 

After  we  have  solved  the  above  modified  auxiliary  parity-condition  problem  (i.e. 
determined   Z  =  Z*   for  (46)  such  that  (34)  holds) ,  we  have  all  the  information  required 
to  determine  without  explicitly  computing  the  entire  force-level  trajectories  whether 
or  not  force  annihilation  occurs  in  battles  modelled  with  (12) .   We  may  apply  Theorem 
3  via  (28)  (possibly  using  (26))  to  see  who  can  be  annihilated  and  use  results  such  as 

(40)  (or  equivalent  ones  for  hyperbolic-like  power  Lanchester  functions  (see  Taylor  and 

r  25 1 

Brown    ))  to  see  if  force  annihilation  occurs  in  a  given  finite   time  (for  example, 

for  battles  modelled  by  (13)).   When  C  =  0   for  the  coefficients  (12)  (e.g.  for  the 
model  (13),   R  =  Ra  s.  R   from  (14)),  we  have  by  Corollary  3.1 

COROLLARY  3.2:   For  combat  between  two  homogeneous  forces  modelled  by  (11) 
and  (12)  with  C  =  0,   the  X  force  can  be  annihilated  if  and  only  if 
Z*xQ  <  (/k  k,  /(y+l))1-2p/k  /k,  yQ,  where  Z*  =  Z*(Y,U,v)   is  such  that 
(34)  holds  for  (46)  . 

We  will  next  give  exact  analytic  results  for  cases  of  no  offset  (i.e. 
A  =  0  =*  y  =  0)   and  discuss  the  difficulties  of  determining  Z*  when  there  is  offset 
(i.e.   A,y  >  0).   Moreover,  results  for  the  special  case  of  no  offset  help  provide 
a  lower  bound  for  Z*  in  the  general  case. 

6.   Results  for  Power  Attrition-Rate  Coefficients  with  No  Offset. 

When  the  offset  parameter  A  =  0,   equations  (46)  become  for  y,v  >   -1 

de  /ds  =  -e^       with  e  (s=0)  =  1, 

—         8  —  — 

de^/ds  =  -s  e      with  e  (s=0)  =  Z. 


(47) 


Solving  (47)  by  successive  approximations,  we  obtain 
e-(s;Z)  -  I      p2ksk(e+2)/{  u  J(j-p)} 

k=0  3=1         CO  k 

-Z  I     p2ksk(e+2)+1/{  *  j(j+P)},  (48) 

k=0  j=l 
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tfhich   is   not   the  most   useful  result    for   large      s.      Without   explicitly   computing     e    (s;Z) 

A 

it  is  impossible  to  determine  from  (48)  how  e  (s;Z)   behaves  for  increasing  s.   However, 

A 

let  us  write   e^CsjZ)   as  [observe  that  e   satisfies  the  generalized  Airy  equation 
(see  SWANSON  and  HEADLEY^190   d2e"  /  ds2  -  sBe~  =  0,  which  is  well  known  to  be  reducible 

A  A 


to  Bessel's  equation] 

e~(s;Z)  =  r(q)p"q{A.(s)+  [l-Zpq_Pr(p)/r(q)][ps1/2I  (S)},  (49) 

A  p  p 

for  any  arbitrary   Z,   where  A  (s)   denotes  the  generalized  Airy  function  of  the  first 

8 

[19] 
kind  of  order  0   (see  Swanson  and  Headley    )  ,   I  (S)   denotes  the  modified  Bessel 

~ —  P 

(8+2^ /2 
function  of  the  first  kind  of  order   p,   and  S  =  2ps       .   Also, 

e^(s;Z)  =  r(q)p"q{(2p/7T)(sinpiT)  s(6+1)/2K  (S)  -  [1-Z  q"Pr(p)  /r(q)  ]dF/ds  (s)  }  ,     (50) 

where  K  (S)   denotes  the  modified  Bessel  function  of  the  third  kind  (also  called 

q 

1/2 
Macdonald's  function)  of  order  q  and  F(s)  =  ps   I  (S) . 

[19] 
The  behavior  of  AQ(s)   for  s  ^  0  is  readily  seen  from  (see  Swanson  and  Headley    , 
p 

A.(s)  =  (2p/Tr)(sinpTT)  s   K  (S)  .   It  is  readily  seen  (see  p.  119  and  p.  123  of  LEBEDEV  *■ 
0  p  

or  pp.  250-251  of  Olver    )  that  K  (S)   is  strictly  decreasing  and  positive,  and 

lim  s  K  (S)  =  0,  where  v   is  any  real  number.   It  is  then  clear  (see  also  p.  1404  of 

reference  19)  that  A  (s)  ,   I  (S)  >  0  V  s  ^  0 ,   lim  A  (s)  =  0,   and   lim  I  (S)  =  +°°. 

p  s->-f°°  s->+°°  P 

Consequently,  the  requirement  that  (34)  holds  (i.e.  that  e   and  ey  behave  like 

strictly  decaying  exponentials)  means  that  the  second  term  in  the  expressions  (49)  and 

(50)  must  vanish.   Thus,  for  y  =  0  we  have 

z*(Y=o,u,v)  =  pp~qr(q)/r(P).  (51) 

Hence,  Theorem  3  becomes  (see  Note  10) 

THEOREM  4:   Consider  combat  between  two  homogeneous  forces  described  by  (11) 
with  attrition-rate  coefficients  (12)  with  A  =  0.   Assume  that  the  model 
applies  for  all  cime  and  that  Y   "wins"  when  x(tf)  =  0  with  y(tf)  >  0. 


Then  Y  wins  if  and  only  if  xQ«>  (s=s  ;Z*)  <  /k  /k  yQ(*/k  k  /(u+1) )  q~Fex(s=s0  ;Z*) 
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where   Z*  =  Z*(Y=0,y,v)  =  pP~qr(q) /r(p)  ,  e~(s;Z*)  =  p~Pr(q)Ae(s)  ,   and 
e~(s;Z*)  =  (2p/7r)(sinpTT)  p_qr(q)s(B+1) '  2K  (S)  .   For  C  =  0,  we  have  s  =  0 
so  that  Y  wins  if  and  only  if  x  r(q)  <  /k  /1c  y  r(p)(/k  k7 /(u+v+2)  )q~P. 

We  observe  that  the  infinite  series  form  (48)  is  not  of  any  value  for  determining 
asymptotic  properties  of  the  solution  e   to  (47)  (and  consequently   Z*) ,  although  it 

A. 

is  useful  for  computational  purposes. 

7.   Offset  Linear  Attrition-Rate  Coefficients. 

When  the  offset  parameter  A  >  0,   explicit  analytic  results  for  Z*     are 
apparently  not  possible.   Before  considering  the  general  case  of   y,v  >  -1  and  y   >   0, 

we  will  find  it  instructive  to  consider  the  special  case  of  offset  linear  attrition- 

m 
rate  coefficients  (i.e.   Y  >  0  with   y  =  v  =  1)  studied  by  Bonder  and  Farrell    .   This 

examination  will  show  us  why  analytic  results  for   Z*  are  elusive  in  cases  of  positive 

offset. 

For  y  >  0  with  y  =  v  =1,   equations  (30)  and  (31)  become 

d2x/ds2-  (I+y/v^x  =  0,  (52) 

with  initial  conditions  x(s=0)  =  x   and  dx/ds(s=0)  =  -/k  /K    yn >  and 

de~/ds  =  -e~                          with  e~(s=0)  =  1, 

de~/ds  =  -(l+Y/^s)e~                    with   e~(s=0)  =  Z. 

Let  us  now  consider  solving  the  above  modified  auxiliary  parity-condition  problem 
(i.e.  determining  Z  =  Z*  for  (53)  such  that  (34)  holds).   Solving  (53)  by  successive 
approximations,  we  may  write 

e^(s)  =  h(s)  -  Zw(s),         and         e~(s)  =  ZH(s)  -  W(s) ,         (54) 

where  h(s)  =  h(s;Y)  =  h(X=s  jU^/v's-)  ,  w(s)  ,   H(s)  ,   and  W(s)   denote  the  auxiliary 

r  25] 
offset  linear  Lanchester  functions  introduced  by  Taylor  and  Brown    .   Infinite  series 

representations  of  these  functions  are  given  in  reference  25.   Subsequent  research  has 

shown  that  these  hyperbolic-like  Lanchester  functions  possess  the  properties  given  in 
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(53) 


Table  I.   Unfortunately,  information  on  the  asymptotic  behavior  of  the  auxiliary  offset 
linear  Lanchester  functions  for  large  s  >  0,   which  is  needed  to  solve  the  modified 
auxiliary  parity-condition  problem  (53)  ,  is  lacking  and  apparently  not  obtainable  by 
the  standard  methods  involving  integral  representation  (see  Ince     or  Olver    ). 
Consequently,  we  have  not  been  able  to  develop  an  explicit  analytic  expression  for 
Z*(Y>0,y=l,v=l) ,   although  we  give  upper  and  lower  bounds  for  Z*  in  the  next  section 
for  general  u,v  >  -1.   Additionally,  we  should  point  out  that  there  are  computational 
difficulties  in  searching  for  Z*  via  its  definition  (34):   (a)  one  doesn't  know  how 
large  to  take  s   for  "satisfactory"  results,  and  (b)  numerical  difficulties  in  evaluat- 
ing e  (s;Z)   and  ey.(s;Z)   as  given  by  (54)  occur  for  large  values  of   s   (since  we 
are  taking  the  difference  of  two  very  large  numbers  and,  at  least  on  a  digital  computer, 
can  retain  only  a  limited  number  of  significant  digits  in  these  numbers) . 

Equation  (52)  is  deceptively  simple  looking.   Using  variation  of  parameters,  we 
may  also  express  its  solution  as 


rs 


x(s)  =  x_  cosh  s  -  y^A.  /k,  sinh  s  +  y 


sinh  (s-a)x(a)//ado .  (55) 


Although  (55)  is  a  simple  looking  expression,  this  Volterra  integral  equation  is,  unfor- 
tunately, no  easier  to  solve  than  (52)  and  leads  to  the  same  results  as  given  by 

[22]  r  251 

Taylor     and  Taylor  and  Brown     .   However,  in  the  next  section  we  show  an  easy 

way  to  obtain  valuable  information  about   Z*  directly  from  (53) . 

8.   Bounds  on  Z*   for  Power  Attrition-Rate  Coefficients  With  Positive  Offset. 

We  will  now  develop  upper  and  lower  bounds  for   Z*(y>0,u,v).   These  bounds 
establish  the  existence  of   Z*   (and  consequently  Q*)   for  general  power  attrition- 
rate  coefficients  (12)  by  the  continuous  dependence  of  solutions  to  (46)  on  the  initial 
conditions . 

The  following  two  lemmas  will  be  used  to  obtain  an  upper  bound  for   Z*(y,y,v) 
for  u,v  >  -1. 
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TABLE  I.   Properties  of  the  Auxiliary  Offset 
Linear  Lanchester  Functions 

1.  dh/ds  =  W,  dW/ds  =  (l+y/^h 

2.  dw/ds  =  H,  dH/ds  =  (1+y//s)w 

3.  h(s)H(s)  -  w(s)  W(s)  =1     Vs 

4.  h(s=0)  =  H(s=0)  =  1 

5.  w(s=0)  =  W(s=0)  =  0 

6.  h(s;y=0)  =  H(s;y=0>  =  cosh  s 
7              w(s;y=0)  =  W(s;y=0)  =  sinh  s 


Note:   We  use  the  notation  h(s)  =  h(s;Y)  =  h(A=s ,u=y//s) 

[251 
(see  Taylor  and  Brown    ),  and  similarly  for  w(s)  ,  H(s) ,   and  W(s) 
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LEMMA   1:      For      6^1      and      x,y   >    0,      26    1{x<S+y5>   ^    (x+y)6. 

r 

PROOF:   For  6^1,   f(x)  =  x   is  a  convex  function.   A  well-known  theorem  for  convex 
functions  says  that   {f (x)  +  f (y) }/2  £  f([x+y]/2),  whence  follows  the  lemma.    Q.E.D. 

LEMMA  2:   For   6*1  and  x,y  >  0,   x5 + y6  £  (x+y)5 . 

PROOF:   Dividing  by   (x+y)  ,  we  need  to  show  that   [x/(x+y)]  +  [y/(x+y)]   £  1.   If  x 
and  y  >  0,   then  x/(x+y)   and   y/(x+y)   are   <  1.   Hence,  for  any  6  £  1  we  have 
[x/(x+y)]   ;>  x/(x+y)   so  that   [x/(x+y)]  +[y/(x+y)]   ^  [x/(x+y) ] +  [y/ (x+y) ]  =  1.    Q.E.D, 

Using  the  above  lemmas,  we  now  prove  Theorem  5.   The  upper  bound  to  be  given 
in  Theorem  5  might  be  improved  upon,  although  we  feel  that  for  computation  determina- 
tion of   Z*  by  interval  search,  it  is  not  essential  to  have  a  better  bound. 

THEOREM  5:   For   u  >  -1  and  Mil,  we  have   Z*(y,y,v)  <  1+2V~  (u+1)2/ 
[(v+l)(u+v+2)]+7V2V_1(u+l)2/(u+2)  ;  while  for  u  >  -1  and  -1  <  v  £  1, 
we  have   Z*(y,u,v)  <  1  +  (u+1)2/  [  (v+1)  (u+v+2)  ]  +  yV(u+l)2/(u+2)  . 

PROOF :   Recalling  (34),  we  have  from  the  first  equation  of  (46)  that  e~(s;Z*)  <  1 

A 

for  s  >  0;   and  from  the  second  equation  of  (46)  we  then  obtain  for  s  >  0  and  v  £  1 

de^/ds  >  -s6(l+Y/sV  *  -2V_1s6(l+YV/sCtV),  (56) 

the  latter  inequality  being  a  consequence  of  Lemma  1.   From  (30)  we  obtain 

e^s)  >  z*-2V-1s(v+1)/(^+1)  •  (vi+l)/(v+l)-2V-V(y+Ds1/(y+1).         (57) 

Using  (57)  and  considering  the  first  equation  of  (46),  we  obtain   e  (s;Z*)  <  U(s;Z*)> 
where 

U(s;Z*)  =  l  +  2v-1{(u+l)2/(u+2)}{s(y+V+2)/(p+1)  •  (u+2)  /  [  (v+1)  (y+v+2)  ] 

+  yVs(m+2)/(u+1)}_z,s> 

Since  we  must  have   0  <  e~(s;Z*)  <  U(s;Z*)   for  all  s  £  0,   it  follows  that  for   s  =  1 

A 

we  must  have  U(s=i;Z*)  >  0,  whence  follows  the  theorem  for   v  £  1.   Lemma  2  and 

similar  arguments  are  used  to  prove  the  theorem  for   -1  <  v  £  1.  Q .E .D. 
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Let  us  now  consider  the  development  of  a  lower  bound  for   Z*(y>0,u,v).   Before 
proving  the  key  lemma  (Lemma  3)  for  the  proof  of  Theorem  6,  we  discuss  some  preliminary 
considerations . 

r  oft  1 

As  shown  by  Taylor  and  Parry     ,  force  annihilation  for  square-law  attrition 
processes  sometimes  may  be  predicted  by  considering  the  force-ratio  equation.   For 
equations  (46),  the  "force  ratio"  u  =  e  /e   satisfies  the  Riccati  equation 

A   i 

du/ds  =  s6(l+r/sa)Vu2- 1,  (58) 

with  initial  condition  u(s=0)  =  1/Z.  We  observe  that  u(sf)  =0   if  X  is  annihilated 
(i.e.   e  (sf)  =  0  but  e^sf)  >  0),   and  u(s-)  =  +°°  if   i"  is  annihilated.   For 
Z  =  Z*  =  Z*(y,y,v)   such  that  (34)  holds,  we  have  a  "draw,"  and  u*(s)  =  u(s;Z=Z*)  >  0 
and  finite  for  all  finite  s  ^  0.   Since  (58)  has  a  singularity  at  s  =  0,  we  apply 

r  •  ,      Ct  ,   s  U+l  -  -  U+l 

the  transformation  t  =  (s  +y)     to  obtain  for  x  ^  Xq  =  y 

du/dt  =  x6u2-  (l-y/Ta)y,  (59) 

with  initial  condition  u(x=y   )  =  1/Z. 

We  now  develop  a  lower  bound  for  Z*(y>0,y,v)   by  comparing  results  for  y   >  0 
with  those  for  y   =   0.   For  y  -   0 ,  we  denote  u  as  w  and  obtain  for  T  ^  0 

dw/dr  =  x6w2-  1.  (60) 


Corresponding  to  Z  =  Z*  =  Z*(y=0,y,v),  we  have  via  (49)  and  (50) 

w*(t)  =  x"B/2K  (T)/K  (T)  >  0,  (61) 

p    q 

(8+2) 1 1 
where  T  =  2px       .   Since  K  (x)   is  finite  and  >  0  for  all  v,x  >  0   (see 
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Lebedev    ,  p.  136),  we  readily  verify  that  w*(t)   is  finite  for  all  T  >  0  and 

that  w*(t=0)  =  pq_Pr(p)/r(q) .   Let  us  observe  that  w(x  )  >  w*(t,)  ■  w(x  ;Z*)  =»  we 

have  w(x)  >  w*(x)  V  x  ^  0  and  w(x_)  =  +00  for  some  finite  i     >   x.  ,   since  D  =  w  -  w* 

o 

satisfies  dD/dx  =  x  (w+w*)D.   Consequently,  w(x;Z)   corresponds  to  Z  <  Z*  and 

w(x;Z)   becomes  infinite  at  some  finite  time  (see  equations  (49)  and  (50)).  We  now 

state  and  prove  the  key  lemma  for  developing  a  lower  bound  for  Z* 
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LEMMA  3:  Let  w*(t)   be  given  by  (61)  and  let  u(t)  =  u(t;Z)  satisfy 
(59)  for  t  ^  t0  =  yy   •   Then  if  u(t.)  £  w*(x  ) ,   it  follows  that 
u(t)  >  0  Vt  ^  :.   and   u(x„)  =  +°°  for  some  finite  x„  >  x..  Conse- 
quently, Z*(y>0,y,v)  >  Z. 

PROOF:   Consider  D  =  u-w.   It  satisfies  for   x  £  xn  via  (59)  and  (60)  the  equation 
dD/dx  =  TB(u+w)D+  {l-(l-Y/xa)y}.   If  D(t1)  =  u(t  )-w(t  )  £  0,   then  dD/dx (t)  >  0 
and  D(t)  >  0  Vt  >  t  ,   Thus,  when  u(x  )  ^  w*(x  )  ,  we  can  find  w(x  )   for  x„  ^  x1 
such  that   u(x„)  >  w(x„)  >  w*(x-)  ,  whence  follows  the  lemma  from  u(x)  >  w(x;Z) 
Vt^  t„  with  Z  <  Z*(y=0,y,v)   and  the  above  observation  that  w(x;Z)   becomes 
infinite.     Q.E.D. 

Letting  x,  =  T_,  we  obtain 

THEOREM  6:   Z*(y>0,u,v)  >  l/w*(Yy+1)  =  t^/2K  (T_) /K  (Tn) ,  where 

U    q   U    p   U 

x„  -  V"+1  a„d  TQ  .  2p,r2,/2- 

Since  w*(x=0)  ^  w*(x)   for  8^0,  we  have  as  an  immediate  corollary 

COROLLARY  6.1:   For  v  £  u,   Z*(y>0,y,v)  >  pP_qr(q) /T(p) . 

Let  us  observe  that  for   3  >  0,   the  lower  bound  given  in  Corollary  6.1  is  weaker  than 
that  in  Theorem  6. 

9.   Future  Computational  Work. 

As  we  have  seen  above  in  Section  4,  force-annihilation  prediction  depends  on 
knowing  the  parity-condition  parameter  Q*,  which  may  be  called  "the  Y   equivalent 
of  an  X  force  of  unit  strength."  We  have  explicitly  determined  Q*   (via  determining 
Z*  for  the  modified  auxiliary  parity-condition  problem  (46))  for  the  power  attrition- 
rate  coefficients  (12)  in  the  case  of  no  offset,  i.e.   A  =  0.   Tabulations  of,  for 
example,  the  new  modified  exponential-like  general  Lanchester  functions   e  (s;Z*)   and 

A 

e^s;Z*)   would  facilitate  force-annihilation  prediction  (see  Theorem  4).   It  remains 
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to  determine  Q*   for  cases  of  positive  offset,  i.e.   A  >  0.  As  discussed  in  Section  7, 
analytic  results  for   Z*  in  the  modified  auxiliary  parity-condition  problem  (46)  with 
Y  =  A*  ( /k  k,  /(u+1))  >  0  are  apparently  not  possible  by  the  usual  analytical 

methods  so  we  must  turn  to  numerical  methods.   It  appears  that  a  large  number  of  cases 
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of  tactical  interest  (see  Taylor  and  Brown  "  )  would  be  covered  by  determining  Z* 

for  u,v  =  0,1,2,3  and  for  a  range  of  values  of  y   >  0.   One  would  be  interested  in, 

for  example,  plotting  Z*  versus  y      for  fixed  values  of  u  and  v. 

Since  we  have  developed  upper  and  lower  bounds  for   Z*  when  y   >  0   (see 

Theorems  5  and  6) ,  we  can  use  standard  one-dimensional  search  techniques  (see ,  for 

r  la] 
example,  WILDE   '  )  to  calculate  an  approximate  value  of  Z*  with  any  predetermined 

degree  of  accuracy  (of  course,  depending  on  how  much  computation  we  wish  to  do).   Since 
(34)  must  hold  for  all  s  ^  0,  we  must  determine  how  long  (i.e.  for  how  large  a  value 
of  s)   to  carry  out  computations  of  e  (s;Z)   and  e^XsjZ)   in  the  modified  auxiliary 
parity-condition  problem  (46)  for  a  given  trial  value  of  Z*   (denoted  as   Z)  to  see 
whether  it  is  too  large  or  too  small.   Although  •  e  <  0  or  e^.  <  0  clearly  terminates 
calculations  for  a  given  Z,   it  would  be  desirable  to  be  able  to  terminate  computations 
before  this  occurs,   especially  when  Z  is  near  Z*.  In  the  future  we  will  show  that 
by  considering  the  Ricatti  equation  (59)  one  can  "cut  off"  conputations  for  a  given 
value  of   Z  well  before  either  of  the  two  annihilation  conditions  is  actually  reached. 
This  acceleration  in  the  termination  of  trial  computations  is  particularly  useful  when 
force  annihilation  occurs  for  large  s   (i.e.  for  Z  near  Z*) .  Additionally,  the 
qualitative  behavior  of  solutions  to  (46)  (or  (11)  ,  for  that  matter)  is  probably  best 
understood  by  considering  the  Ricatti  equation  (59)  (or,  equivalently ,  (58)). 

As  discussed  at  the  end  of  Section  4,  prediction  of  force  annihilation  within 
a  given  finite  time  involves  the  use  of  tabulations  of  the  quotient  of  two  linearly- 
independent  general  Lanchester  functions.  We  have  indicated  in  Section  4  that  the 
hyperbolic-tangent-like  Lanchester  functions  (e.g.   T  (t)   as  defined  by  (42))  are  to 

A. 

be  preferred  for  reasons  of  the  accuracy  of  their  numerical  computation.   Thus,  there 
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is  a  need  for  tabulations  of  T  (t)  ,  whose  range  is   [0,1/Q*)   for   te  [0,+°°).   For 

the  power  attrition-rate  coefficients  (12)  with  no  offset,  the  power  Lanchester  functions 

(also  called  Lanchester-Clif f ord-Schlaf li  (or  LCS)  functions),  however,  were  inappro- 
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priately  defined  in  Taylor  and  Brown     to  yield  such  tabulations.   Thus,  our  newer 

theory  of  force-annihilation  prediction,  which  also  involves  tabulations  of  canonical 

solutions  (i.e.  canonical  Lanchester  functions)  to  variable-coefficient  Lanchester-type 

equations  of  modern  warfare,  has  suggested  some  refinements  in  the  definition  of  Taylor 
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and  Brown's     auxiliary  power  Lanchester  functions  (see  Note  11).   It  would  be  desir- 
able then  to  redefine  the  LCS  functions  to  fit  within  the  framework  of  Section  A  and 
to  develop  tabulations  of  their  quotient.   If  this  were  to  be  done,  linear  combat  models 
with  power  attrition-rate  coefficients  (no  offset)  could  be  analyzed  with  somewhat  the 
same  ease  as  constant-coefficient  linear  models  (i.e.  (1)). 

10.   Extension  to  a  More  General  Model. 

In  this  section  we  show  that  all  our  above  force-annihilation  results  (except 
those  for  force-annihilation  within  a  fixed  finite  time)  also  apply  to  a  special  case 
of  a  more  general  model.   Moreover,  comparison  techniques  may  be  used  to  extend  these 
results  in  weakened  form  to  the  general  case  of  this  more  comprehensive  model. 

Let  us  consider  the  following  Lanchester-type  equations 


dx/dt  =  -a(t)y-g(t)x  with  x(t=0)  =  x  , 

dy/dt  =  -b(t)x-a(t)y  with  y(t=0)  =  yQ , 


(62) 


where  a(t)  ,   b(t),   ct(t)  ,   and  g(t)  >  0.   We  may  think  of  these  equations  as  modelling, 
for  example,  aimed-fire  combat  between  two  homogeneous  infantry  forces  with  superimposed 

effects  of  supporting  weapons  (which  are  not  subject  to  attrition  and  deliver  area 

r  jf.  I 
fire  against  the  enemy  infantry)  (see  Taylor  and  Parry     for  a  further  discussion 

of  this  model  (62)).   In  this  case,   a(t)   and  8(t)   are  attrition-rate  coefficients 
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which  reflect  the  fire  effectiveness  of  the  supporting  weapons     .   Then,  the  force 

ratio  u  =  x/y  satisfies  the  generalized  Riccati  equation 
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du/dt  -  b(t)u2+  {o(t)  -  B(t)}u-a(t)     with  u(t=0)  =  xQ/y  .  (63) 

For  equal  effectivenesses  of  the  supporting  fires   [i.e.   a(t)  =  0(t)],   equation 
(63)  simplifies  to 

du/dt  =  b(t)u2-a(t),  (64) 

which  is  the  same  Riccati  equation  satisifed  by  the  force  ratio  for  the  model  (11). 
Hence,  when  a(t)  =  3(t)  V  t  iO,   a  battle's  outcome  (in  terms  of  the  force  ratio) 
is  the  same  for  the  two  models  (11)  and  (62)  ,  although  the  battle  ends  more  quickly 
for  (62).   Thus,  in  this  special  case  all  our  above  results  on  force  annihilation  with- 
out time  limitation  (e.g.  Theorems  2  through  4  and  their  corollaries)  developed  for  (11) 
[in  general  or  with  the  coefficients  (12) ]  also  apply  to  the  more  general  model  (62) . 
Furthermore,  comparison  techniques  (see,  for  example,  Hille    )  may  be  used  to  extend 
these  results  in  weakened  form  to  (62).   Consequently,  we  see  that  the  force-annihi- 
lation results  developed  in  this  paper  are  indeed  of  a  fundamental  nature. 

11.   Summary. 

We  have  presented  a  mathematical  theory  for  predicting  force  annihilation  for 
variable-coefficient  Lanchester-type  equations  of  "modern  warfare"  for  combat  between 

two  homogeneous  forces  without  explicitly  computing  force-level  trajectories  (see 

[25] 
Note  12).   This  theory  both  generalizes  and  complements  Taylor  and  Brown's     theory 

of  canonical  solution  forms.   Our  force-annihilation  theory  provides  guidance  for 

certain  parameter  determinations  and  development  of  tabulations  of  Lanchester  functions 

(beyond  those  suggested  in  reference  25)  that  would  allow  one  to  parametrically  analyze 

variable-coefficient  models  with  somewhat  the  same  facility  as  constant-coefficient 

ones. 

We  have  shown  that  force  annihilation  can  be  predicted  from  initial  conditions, 

without  explicitly  computing  force-level  trajectories,  by  knowing  a  parity-condition 

parameter  Q*,   which  is  the  solution  to  a  canonical  auxiliary  parity-condition  problem. 

In  general,  this  prediction  would  be  facilitated  by  having  tabulations  of  certain 
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Lanchester  functions  available.   In  particular,  to  apply  Theorem  3  it  would  be  convenient 
to  have  tabulations  of  the  new  exponential-like  general  Lanchester  functions  E   and 

A 

eT.   However,  if  one  wishes  to  invoke  the  equivalent  Theorem  3',  then  tabulations  of 
the  hyperbolic-like  general  Lanchester  functions  C  ,   S  ,   C  ,   and  S   would  be 

A     A     I  i 

desirable  to  have.   The  parity-condition  parameter  Q*  was  shown  to  be  related  to 

the  range  of  the  quotient  of  two  hyperbolic-like  general  Lanchester  functions  introduced 

[251 
by  Taylor  and  Brown    .   Consequently,  our  force-annihilation  theory  not  only  provides 

new  information  about  the  mathematical  properties  of  hyperbolic-like  Lanchester  functions 

but  also  provides  guidance  for  selecting  canonical  Lanchester  functions. 

We  applied  our  general  theory  to  the  specific  case  of  general  power  attrition- 
rate  coefficients.   Considering  a  modified  auxiliary  parity-condition  problem,  we 
explicitly  determined  Q*   (via  Z*  of  the  modified  problem)  for  power  attrition-rate 
coefficients  with  no  offset  and  gave  upper  and  lower  bounds  for   Z*   for  cases  of 
positive  offset.   Consequently,  in  the  future  one-dimensional  search  techniques  may 
be  used  to  numerically  determine  approximate  values  of   Z*  as  a  function  of  the  off- 
set parameter  A.   We  finally  showed  that  certain  of  our  force-annihilation  results 
also  applied  to  a  more  general  linear  differential  equation  combat  model. 

These  results  may  be  used  in  the  analysis  of  the  dynamic  combat  interactions 
between  two  homogeneous  forces  with  time-  (or  range-)  dependent  weapon  system  capa- 
bilities.  There  is  interest  today  in  such  analytic  models  because  of  improved  opera- 
tions research  techniques  for  predicting  Lanchester  attrition-rate  coefficients,  in 

particular  their  temporal  variations  (see  references  2  through  5  and  25) .   Further 

[41 
discussion  of  such  applications  may  be  found  in  Bonder  and  Farrell    ,  Bonder  and 
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Honig1  \    Taylor1    ,  and  Taylor  and  Brown 
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NOTES 

1.  Bonder  and  Honig    point  out,  however,  that  force  annihilation  may  not  be  the 

appropriate  criterion  for  evaluating  many  military  operations,  expecially  when  force 

[4] 
annihilation  does  not  occur.   See  pp.  192-242  of  Bonder  and  Farrell    for  a  detailed 

Lanchester-type  analysis  of  an  attack  scenario  for  which  other  "end  of  battle  condi- 
tions" play  the  major  role  in  the  evaluation  process.   Nevertheless,  it  is  of  interest 
(especially  for  developing  insights  into  the  dynamics  of  combat)  to  be  able  to  easily 
predict  (without  computing  force-level  trajectories)  the  occurrence  of  force  annihilation 
Such  results  are  not  only  of  intrinsic  interest  but  also  are  useful  in  the  optimization 

of  combat  dynamics,  i.e.  the  combining  of  Lanchester-type  models  with  generalized 
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control  theory  (i.e.  optimal  control/differential  game  theory) (see  Taylor         ). 

2.  Previous  work  by  Bonder  and  Farrell    ,  Taylor     ,   and  Taylor  and  Brown     shows 
that  new  transcendental  functions  arise  even  in  the  case  of  linear  attrition-rate 
coefficients  reflecting  weapon  systems  with  different  effective  ranges  (i.e.  the 

coefficients  (12)  with  u  =  v  =  1  and  A  >  0) .   For  example,  the  differential  equation 

[12 
(52)  could  not  be  found  among  the  445  linear  second  order  equations  tabulated  in  Kamke 

Moreover-  even  when  one  can  express  a  solution  in  terms  of  previously  known  transcendents 
the  appropriate  tabulations  (see,  for  example,  ABRAMOWITZ  and  STEGUN   )  may  not  exist 
(see  Section  5  of  reference  25) .   As  the  equations  of  mathematical  physics  have  pro- 
vided interest  in  many  previously  studied  transcendents,  so  does  the  Lanchester  theory 
of  combat  with  time-varying  fire  effectivenesses  provide  interest  in  new  transcendents. 

3.  The  influential  19th-century  German  military  philosopher,  Carl  von  Clausewitz 
(1780-1831),  stated  in  his  classic  work  On  War  (Vom  Kriege)  (see  p.  276  of  reference  7), 
"The  best  Strategy  is  always  to  be  very  strong,  first  generally  then  at  the  decisive 
point.   ...There  is  no  more  imperative  and  no  simpler  law  for  Strategy  than  to  keep  the 
forces  concentrated. " 
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4.  As  pointed  out  in  reference  26,  the  entire  topic  of  modelling  battle  termination 
is  a  problem  area  in  contemporary  defense  planning  studies,  and  there  is  far  from 

universal  agreement  as  to  even  which  variables  should  be  taken  as  the  significant 
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variables  for  modelling  this  complex  process.   For  further  references  see  Taylor 

[4] 

5.  Bonder   and   Farrell  take   range    (i.e.    force   separation)    to  be   the   independent 
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variable   in   their  work,   while  Taylor  and   Taylor   and   Brown  take   time  as  we  have 

done  in  this  paper. 

6.  Recalling   the   constant   coefficient   result    (6) ,   we   consider 

d{exp(/ab   t)  }/dt   =  a/b/a  exp(/ab   t)    =  b/a/b   exp(/ab   t)  , 
and 

d{exp(-v/ab   t)  }/dt   =   -a/b/a  exp(-/ab~  t)    =   -b/a/b   expC-i/ab"  t)  , 

to  obtain  motivation  for  (15)  and  (16) . 

7.  We  take  a(t)   and  b(t)   to  be  analytic  in. the  (finite)  complex  plane  except  for 
a  finite  number  of  singularities  on  the  real  axis.   The  singularities  of  (17)  then 
occur  at  the  zeros  and  singularities  of  a(t)   and  at  the  singularities  of  a(t)b(t). 

Consequently,   t_  belongs  to  the  set  of  points  consisting  of  the  zeros  and  singulari- 
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ties  of   a(t)   and  b(t)  (see  Taylor  and  Brown    ).   We  define   t    this  way  in  order 

to  reduce  the  number  of  tabulations  of  exponential  general  Lanchester  functions  required 
for  force-annihilation  analyses  (see  Theorem  3).   For  example,  for  the  general  power 
attrition-rate  coefficients  (12)  we  have   t   =  -C ,   and  for  fixed  A  £  0 ,   u,   and  v 
only  a  single  tabulation  of   e    and  e^      is  required  to  handle  all  problems  with 
C  ^  0   (see,  for  example,  Theorem  4). 

8.  It  is  well  known  (see,  for  example,  pp.  647-650  of  HilleL    or  p.  120  of  Kamke 
that  the  quotient  of  two  linearly  independent  solutions  to  (31)  ,  which  is  equivalent 

ri8i 

to  (17),  satisfies  Schwarz's  (third  order)  differential  equation  (see  SCHWARZ    ) 
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Cn,s}  =  -2i(s), 

where  n  denotes  the  quotient  of  two  linearly  independent  solutions  to  (31)  [i.e. 
n.  =  ev(s)/ev(s)  ] ,  {r),s}  =  n1"  -  (3/2)  (n")  2/n '  denotes  the  Schwarzian  derivative  of 
n.  with  respect  to  s,  n1  denotes  dn/ds,  etc.,  and  I(s)  denotes  the  invariant 
of  the  normal  form  (31).  For  numerical  computation  of  n. ,  however,  there  is  no  advan- 
tage to  consider  this  third  order  equation,  and  it  is  preferable  to  calculate  n  from, 
for  example,   n(s)  =  e  (s)/e  (s),  where  e   and  e   also  satisfy  (29)  and  (30). 

A       A  A  A 

9.  Recalling  our  development  of  (7),  we  see  that, except  for  the  special  case  in  which 

(23)  holds,  this  same  approach  fails  to  yield  the  time  for  X  to  reach  his  breakpoint 

BP  BP 

(assumed  to  be  positive)  [i.e.   t    such  that  x(t=t   )  =  x^p  >  0] .   Consequently,  it 

is  apparently  impossible  to  predict  in  the  manner  described  in  the  main  text  the  out- 
come of  a  fixed  force-level  breakpoint  battle  with  positive  breakpoints  unless  (23) 
holds. 

10.  For  the  case  of  power  attrition-rate  coefficients  with  no  offset  (i.e.   A  =  0   in 
(12)),  the  second  annihilation  condition  given  in  Theorem  4  (i.e.  the  one  for  C  =  0) 
and  an  equivalent  form  of  the  first  (i.e.  the  one  for  C  ^  0)  may  be  developed  by 
inspection  when  one  expresses,  for  example,  the  time  history  of  the  X  force  level 

(which  satisfies  (17))  in  terms  of  the  so-called  generalized  Airy  functions  (see  Swanson 

r  19 1 

and  Headley    ).   For  example,  for  C  =  0  we  have 

x(t)  =  (pP/(2p)){X+Y}A8(T)  +  (pP/(2/p~)){X-Y}BB(T), 

where  A„   and  B„  denote  the  generalized  Airy  functions  of  the  first  and  second  kinds 
of  order   g,   X  =  xQr(q),   Y  =  YQ^J\    r(p)  (/k^/ (u+v+2) ) 1_2p ,   and   T  = 
(A  K  /(u+1))  Pty   .   The  result  given  in  Theorem  4  for  C  =  0  follows  from  the  proper- 
ties of  the  generalized  Airy  functions  (i.e.   A  (O ,  B  (?)  >  0  V  |  ^  0,   lim  A  (?)  =  0, 

v      v  ?-H<» 

and  lim  B  (?)  =  +»)   and  the  above  representation  for  x(t) .   Unfortunately,  this  result 

does  not  generalize  to  other  cases  of  interest,  although  it  did  motivate  our  general 
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theory  of  force  annihilation  developed  in  this  paper.   The  generalized  Airy  functions 
may  be  considered  to  be  generalizations  of  the  exponential  function  (see  p.  446  of 
reference  1  or  p.  393  of  reference  17  for  plots  of  the  standard  (i.e.   8=1)  Airy 
functions)  and  arise  in  the  study  of  the  asymptotic  behavior  of  solutions  to  certain 
differential  equations. 

11.  Although  theoretically  results  for  power  attrition-rate  coefficients  with  no 

offset  are  expressible  in  terms  of  "known"  transcendental  functions,  new  Lanchester 

[251 
functions  were  introduced  by  Taylor  and  Brown     because  of  lack  of  tabulations  of 

these  in  many  cases  of  interest. 

[91 

12.  In  his  well-known  survey  paper  on  the  Lanchester  theory  of  combat,  Dolansky 

suggested  the  development  of  outcome-predicting  relations  without  solving  in  detail 
as  one  of  several  problems  for  future  research. 
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